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relativistic limit. A necessary and sufficient condition for a black hole limit is, that 
the electric potential in the co-rotating frame is constant on the disc. In that case 
certain other relations follow. These relations are reviewed with a highly accurate 
post-Newtonian expansion. 

Remarkably it is possible to survey the leading order behaviour close to the black 
hole limit with the post-Newtonian expansion. We find that the disc solution close to 
that limit can be approximated very well by a “hyperextreme” Kerr-Newman solution 
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1. Introduction 

Bardeen and Wagoner [1] were able to solve the general-relativistic problem of a 
uniformly rotating disc of dust by means of a post-Newtonian expansion to very high 
order. Although their method involved the successive numerical evaluation of integrals 
they achieved an impressive accuracy even in the extreme relativistic limit. From the 
exterior perspective, the spacetime approaches the metric of an extreme Kerr black hole 
(outside the horizon) in that limit. This remarkable prediction was rigorously confirmed 
with the exact solution to the disc problem by Neugebauer and Meinel |i2j, see also pj. 
Starting from the exact solution, Petroff and Meinel [4] presented an analytic post- 
Newtonian expansion scheme up to arbitrary order. Recently |5], a similar analytic 
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expansion scheme was found for the more general problem of a uniformly rotating disc 
of electrically charged dust (with a constant specific charge) - without knowing the 
corresponding exact solution. Again, evidence for a black hole limit was given. Based 
on a further elaboration of the method presented in [5], the aim of the present paper is 
a detailed investigation of the electrically charged disc of dust in the extreme relativistic 
limit. In particular, the behaviour of the gravitational and electromagnetic multipole 
moments is studied. The results are fully consistent with a spacetime that approaches, 
from the exterior perspective, the metric of an extreme Kerr-Newman black hole in this 
limit. 

The paper is organised as follows. Section [2] is an introduction to rotating charged 
dust and its mathematical description. Here we also give a proof for an important 
parameter equation. In section [3] we give an outline to the post-Newtonian expansion 
of rigidly rotating discs with constant specific charge and how we can treat the question 
of its convergence. Section [4] is dedicated to the black hole limit of these discs, and in 
section [5] we make statements on the leading order behaviour close to that limit. In the 
Appendix we have provided the first few analytic coefficients of the expansion for the 
relevant quantities within this paper. 


2. Description of rotating configurations of charged dust in equilibrium 

2.1. Basic field equations and model of matter 
The Einstcin-Maxwell equations reac(|] 

Rab ~ ^Rg a b = 87r T ab , (1) 

F[ab;c] = 0 , ( 2 ) 

F ab ,b = 4tt/. (3) 

The homogeneous Maxwell equations (|2]) are fulfilled by a four-potential A a with 

F a b = A b . a — A a . b = A ba — A ab . (4) 

The energy momentum tensor T ab for the model of charged dust that we use can be 
expressed as the sum of a dust part and an electromagnetic part 

T* = + l£“> = ^u b + 4- (f^ - \ gab F cd FA . (5) 

The baryonic mass density p is related to the charge density p e i via 

Pei = ep, (6) 

where in general the specific charge e is a free function. We assume a purely convective 
four-current density, i.e. 

f = p e \u a = epu a . (7) 


i We use geometrized Gauss units c = G = 47reo = 1 and the metric signature (+, +, +, —). Covariant 
differentiation is denoted by a semicolon and partial differentiation is denoted by a comma. 
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Now we consider isolated and rotating configurations of charged dust in equilibrium, 
which leads to a stationary spacetime, giving a Killing vector Furthermore, assuming 
axisymmetry, we have a second Killing vector 77 that commutes with For an 
asymptotically flat spacetime, the far field behaviour of the Killing vectors is given 
by 


re 


-1 + 


2 M 
r 


+ O (r- 2 ) , 


f a Va 

rf'Tjh 



+ O (r- 4 ) , 


( 8 ) 


with r = \Jx 2 + y 2 + z 2 and x, y, z being asymptotic Cartesian coordinates, the 
gravitational mass M and the angular momentum J. The metric can be written globally 
in terms of Lewis-Papapetrou-Weyl coordinates 

ds 2 = f^ 1 [h (dp 2 + d£ 2 ) + p 2 d(p 2 ] — / (dt + adip) 2 , (9) 


with 


e = d tl 77 = d v . (10) 

For a more detailed discussion of equilibrium states and the far field see |6j. Using 
Lorenz gauge, we get A a = (0, 0, A^, A t ). 

Altogether, we have to solve the Einstein-Maxwell equations for the five unknown 
functions /, h, a, A t and A^ which depend only on the coordinates g and f while the 
angular velocity U = ^ and the specific charge e are specified by a particular physical 
model. 

Here we choose rigid rotation (f2 = constant) and a constant specific charge e. In 
our units we have e G [—1,1], obtaining static electrically counterpoised dust (ECD) con¬ 
figurations (see [7] and ii) for e = ±1. In the case of rigid rotation, a global co-rotating 


frame exists with ip' = cp — Qt, where the metric retains its form 

ds /2 = f'~ l \h' (dp /2 + d<C /2 ) + p /2 dy/ 2 ] — f (dt' + a'dy/) 2 (11) 

and the metric functions and electromagnetic potentials transform as 

/'= / (1 + fla) 2 - -j-, (1 + fla) / = (1 — Qa) /', — = — (12) 

A' t , = A t + flA v , A!^, = A<p. ( 13 ) 

The four velocity in the co-rotating frame can be expressed as 

u' n = ^0, 0, 0,1/ v / 7') i <= (o.O.-a'^-A/T 7 )- ( 14 ) 


2.2. Ernst equations 

Outside the dust configuration, we apply the Ernst equations [9] which are equivalent 
to the stationary and axisymmetry electro-vacuum Einstein-Maxwell equations. Using 
the abbreviation 


A t = -a, 


(15) 
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they reac^§] 

(fR£ + <5$) AS = (VS + 2$V4>) • VS, 
(RS + <&$) A4> = (VS + 2i>V4>) ■ V<h 


with 


and 


S = /-$$ + \b, 4> = a + i/3 


(16) 

(17) 

(18) 


£2 £2 

6,e = a,c - 2 «Af> + 2£a e , b c = — a e - 2a/3 jC + 2£a c , (19) 

£? 


/ / 

Ae = — ( aa ,C d* ^A,c) 5 Ac = -( ac he + Ay )Q ) . 


( 20 ) 


Here the differential equations for the metric function h are decoupled from the others. 
It can be calculated via a path-independent line integral afterwards: 


( ln Aa = \e (W)i - ( ln /)* c - ^ - 4)) + - <c) 

~ 2 ~—^ {At, e - - 4-^ (H^H^ - , 

(Inh) ^ = Q ^(ln/) e (ln/) iC - ^a iQ a^ + A^A VtB A ViC - 4 g ^ A tjB A tx 


( 21 ) 


r 

— 4— + A^H^). 

Q 


( 22 ) 


A 5. Multipole moments 

We introduce the potentials 
l-S 


and q = 


2<f> 


1 + £ 1 + S 

On the positive axis of symmetry they can be written as a series expansion in ( 


(23) 


- 1 . 


OO 

1 ^-\ m r 

C ( n 

s n =0 S 


OO 

1 e r 

^ = 72^0 

s n =0 S 


(24) 


The gravitational multipole moments P n and the electro-magnetic multipole moments 
Q n introduced by Simon mu can be obtained from the coefficients m n and e n (see HU 
and nzn- For n < 3 they simply read P n = m n and Q n = e n . 


2 . 4 . Physical quantities 

In [T2j we presented the parameter equation 

M = 20 J + VM 0 (25) 


§ We use the operators V and A as in Euclidean 3-space, with cylindrical coordinates (g, (, ip). 
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for the gravitational mass. Here we will give a short derivation. The Lorentz force 
density is given by 

f“ = - (T (em)ab ). 6 = F ab j b (26) 

and orthogonal to the Killing vectors 

- £ a f “ = (£ a T^ ab ) ]b = 0, - Va f a = ( Va T^ ab ) ]b = 0. (27) 

Gauss law, the local conservation of energy momentum and (27) give us the option to 
define a set of physical quantities in a coordinate independent way. 

Gravitational mass M and angular momentum J are given by (see also El) 


M = 2 


(Uj + i iigfj Cn‘dV, J = ~j WnMV 


(28) 


V V 

where V is a space-like hyper-surface with the 3-dimensional volume element dV = 
y/Oigd 3 x and the future pointing unit normal n a . The definitions of the gravitational 
mass M and the angular momentum J are arranged in such a way, that they coincide 
with the far held behaviour in (|8j). From the local mass conservation law (y [ ru a ). a = 0 
and the continuity equation f. a = 0, we get the baryonic mass M 0 and the total charge 
Q as 


Mn — / dM n — 


Vdust 


[iu a n a dV, Q = d Q = 


efiu a n a dV. (29) 


Vdust 


Vdust 


Vdust 


With (27) we can define the electromagnetic held energy M em and the electromagnetic 
held angular momentum J em as 


M em = 


vywiv. - - / r^r^dv. 


(30) 


v 


V 


The dehnitions in (30) are arranged in a way, that they coincide with the electromagnetic 
held energy and the electromagnetic held angular momentum in classical electro¬ 
dynamics. 

In Lewis-Papapetrou-Weyl coordinates we integrate (30) by parts and use the 
inhomogeneous Maxwell equations ([3]). This leads to 


M em = - 


A aJb (C-k} v a )n b dV, J em = - / A ajb rj a n b dV. (31) 


^dust 


^dust 


So the domain of integration V in (30) (and therefore also in (28)) can be reduced to 
the volume of the dust. For rigid rotation we can combine (31) and (28) with the help 

M — 2QJ = [ (y/f' + ea') dM 0 


of (flil) to 


(32) 


Vdust 


where f and the electric potential of = —A' t , refer to the co-rotating system. For a 
constant specihc charge, it follows from the equations of motion 

/ih a = f a , 


(33) 
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evaluated in the co-rotating frame, that 
\[Y' + col — V — constant 
inside the dust and so we get 

M = 2 QJ + VM 0 

independent from the shape of the dust configuration. 


(34) 

(35) 


3. Post Newtonian expansion of rigidly rotating discs with constant 
specific charge 

3.1. Mathematical problem 

The dust configuration shall be a circular disc with the coordinate radius go, centred 
in the origin of the Lewis-Papapetrou-Weyl coordinate system (£ = 0, 0 < g < go). 
The mass density /i can be written formally as p = a p (g)\/ f /h5(() with the proper 
surface mass density < 7 p (p). By integrating the full Einstein-Maxwell equations over 
a small flat cylinder around a mass element of the disc, matching conditions between 
the metric functions and electromagnetic potentials (and their derivatives) at f = 0 ^ 
(above and below the disc) can be derived. The interior geometry of the disc itself 
and the surface mass density can be calculated from these functions as well. Assuming 
reflectional symmetry, we showed in [5] that the whole problem can thus be reformulated 
as a boundary value problem to the Ernst equations (16) and ( |I7| ) with the boundary 
conditions 

= 0, 3$' = 0, £ 0 (W + e y 7 / 7 ) . = 0, £ 0 (v 7 / 7 + ? = 0 (36) 

on the disc. The last one is equivalent to (34). We also presented a powerful post- 
Newtonian expansion in the parameter 

7 = 1 - VX wh ere f c = f{g = 0, C = 0) (37) 

and an algorithm to solve the equations to arbitrary order in 7 analytically as well 
as some results for the Erst eight orders. Here we extended the calculations to the 
tenth order, where the general structure of the solution (global pre-factors, polynomial 
structure, etc.), described in [5], is still valid. 


3.2. New charge parameter 

The disc solution is a three parameter solution L(/j s ,pi ,p 2 ), with p s as a scaling 
parameter. Here we will use the coordinate radius of the disc po or the gravitational 
mass M as scale parameters. Functions and coordinates scaled with g 0 are labelled with 
a * and those scaled with M are labelled with a °. The other two parameters are 7 
(strictly speaking ^y) and e = Q/M 0 . For the purpose detailed in section |4| it is more 
effective to use the dimensionless parameter if = Q/M instead of e, that is given as 


■0(7, e) = e { 1 + i 1 - e2 ) 


0 + ik< 21 


ll£ 2 ) 7 2 + o ( 7 3 ) 


(38) 
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containing e as a global pre-factor. In the ECD-limit (e = ±1) we have M = Mq and 
therefore if G [—1,1]. Functions with a global pre-factor \/l — e 2 also have a global 
pre-factor a /1 — " 0 2 - 


3.3. Question of convergence 


A clear statement to the radius of convergence for a power series (39) could only be 


made, if all coefficients up to infinity, or the full analytic expression are known. 


F{ 7,e) = -Pnl(T) e ) 


OO 

l + ^2cn(e)f 

n=l 


(39) 


with F NL ( 7 ,e) describing the Newtonian limit. Real calculations could only be done to 
a finite value of n: 


F{j,e,N) = F nl ( 7 ,e) 


N 

1 + ^2c n (e)Y' 

n =1 


(40) 


Here we know the first nine coefficients (N = 9). For N > 4 all terms in the expansion 
of the Ernst equations (16) and ( pA| contribute to the result. So it seems very likely, 
that the first nine coefficients c n (e) indicate the general behaviour and in particular the 
convergence of the expansion functions from the post-Newtonian expansion. 

Nevertheless, we can review the situation from a physical point of view. The 
expansion parameter is strictly speaking (see the equations (34) in [5]). So we want 
to know if we have a well defined physical system in the parameter space yjf G [—1,1] 
and e G [—1,1]. All functions are either even or odd in ^7 and either even or odd 
in e. Functions related to the rotation of the disc are odd functions in and have 
the global prefactor \J\ — e 2 . (Note that = ^\J 1 — e 2 .) Functions related to the 
charge of the disc are odd functions in e. In the limit —> +1 we will approach the 
ultra-relativistic limit. Changing the sign of ^7 to negative values results in a disc 
rotating in the opposite direction. As the potentials b and [3 will change their sign, we 
get the complex conjugated Ernst potentials that satisfy the Ernst equations as well. 

We can also describe the parameter transition from positive to negative values of 
y/j. For small values of y'y we reach the Newtonian limit. I 11 the limit y'y —>- 0 the 
angular velocity and the mass density /a will go to zero. The global solution is the 
Minkowski spacetime. For small negative values of ^7 we reach the Newtonian limit 
of the disc with a negative angular velocity D. So we have a smooth transition from 
positive to negative values in Ak /y. 

While the sign of y'y is related to the direction of the rotation of the disc, the 
sign of e is related to the sign of its charge. For e = ±1 we reach the static ECD- 


limit and the equations of motion (33), which can be interpreted as a generalised force 


balance for every dust particle inside the disc, can now be interpreted as the balance of 
electric repulsion and gravitation (see ra)- For | e | > 1 the electric repulsion cannot be 
compensated, so the spacetime would not be stationary anymore. 
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In summary, from the physical point of view it is possible for the post-Newtonian 
expansion to converge for all Ak /y < 1 if |e| < 1. It is sufficient to inspect only positive 
values of v /y and 


e. 


A sufficient condition for the convergence of the power series (39)for 7 G [0,1] is the 

OO 

generalised Raabe-Duhamel test: An infinite series c n (e) is absolutely convergent if 


n =0 


= n 1 


A1+1 (^) 


> 1 , ee[0,l] 


(41) 


for all n > jiq. It was used by Bardeen and Wagoner in [1] describing the uncharged 
case. There the c n are simply numbers, while here we have to deal with functions in 
e or if. So (41) is useful in many cases, but it fails if the c n (e ) have zeros in e which 
depend on n. This does not necessarily mean that the power series is not converging, 
cf. Figure |4j 

Here we can check the convergence criterion (41) only up to n = N = 9. So it can 
be used to support the assumption that N = 9 > n Q for the most quantities from the 
post-Newtonian expansion we have studied here. This is because the K n are increasing 
with increasing n. In many cases the increasing of the IK n is most slowly for e = 1, 
which is the limit of stationarity. 

As a last support for the convergence of (40) we investigate its Pade approximation, 
which gives highly accurate results, as compared with the exact solution [ 2 ] for the 
uncharged case even in the ultra-relativistic limit, see [4.]. 


3-4- Pade approximation 


In many cases the convergence of a power series could be accelerated by using the Pade 


approximation (see USD- The Pade approximation to ( |40| ) is given by 

(e)yV 


P [L, 2JV - L] {F( 7 , e, JV)} = F NL ( 7 , e) 


2 N-L 

1+ E Pn 

n= 1 


1 + X) S n( e )\/7 r 

72—1 


(42) 


Therein, the coefficient functions p n (e) and s n (e) could be determined by the c n (e). For 


that purpose we can equate the coefficients in from the series expansion of (42) with 


(40). For the Pade approximation we use y'y to have a wider scope in choosing L to 


avoid or minimize the numbers of zeros in the denominator of (42) 


4. The black hole limit 

Here we will collect results from the post-Newtonian expansion that indicate most 
strongly the transition to an extreme Kerr-Newman black hole at 7 —> 1. 

It was shown in iCj that a necessary and sufficient condition for a black hole 
limit of a uniformly rotating (uncharged) fluid body is given by M — 2 flJ —>■ 0. In such 
a limit, the square of the Killing vector x = £ + tends to zero on the surface of 
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the body: xfXi 0- From the exterior perspective, this corresponds to the defining 
condition of a stationary and axisymmetric black hole. On its horizon hi, the Killing 
vector x = + 0?7 is null, with O being the angular velocity of the horizon. In the 


uncharged disc case (e = 0), this can immediately be seen from (34) and (35). Note 


that f = —yfXi- With charge (e ^ 0) a black hole limit of our rotating disc occurs if 
and only if, for 7 —y 1 , the electric potential in the co-rotating frame of = — xfAi tends 
to a constant on the disc: 


7 —y 1 : ol(q* < 1 , C* — 0 ) = constant, 


(43) 


since 7 = 1 and of = constant together with (34) mean f = constant = /' = f c = 0 on 

and the relation 


the disc. Because of the Kerr-Newman black hole uniqueness, see 
M = 2 O J + a'Q 


(44) 


following from (35), we are inevitably led to an extreme Kerr-Newman black hole, 
cf. 


There is also another simple argument showing that the black hole limit of the 
disc must lead to an extreme black hole: For the Kerr-Newman black hole, the horizon 
hi is located at a finite interval on the £°-axis {(f — 0 , |£°| < /), that turns to zero in 
the extreme case (1 = 0). This is the only possible “compromise” with the shape of the 
disc, and it means in addition, that the coordinate radius £> 0 of the disc must tend to 
zero in the limit: 


7 —^ 1 : £0 0. 


(45) 


It should be noted that we are dealing here with the exterior perspective, where the 
disc shrinks to the origin of the Weyl coordinate system, which also gives the location 
of the Kerr-Newman black hole’s horizon. As discussed by Bardeen and Wagoner (for 
the uncharged case), the disc region can be considered as a “singularity” in the horizon, 
whereas from the interior perspective (corresponding to finite coordinates g*, (*), the 
disc remains regular for 7 —>■ 1 and is surrounded “by its own infinite, nonasymptotically 
flat universe” p. For a detailed discusion of these two perspectives in the ECD case, 
see [ 8 ]. 

In coordinates normalized by the gravitational mass M the extreme Kerr-Newman 
black hole can be expressed with the parameter if. The Ernst potentials on the axis of 
symmetry are given as 


_ C°- 1-iy/l^W 
C° + 1 — i\/l — if 2 ' 

= _ ^ _ 

C° + 1 — i^/l - '0 2 


(46) 

(47) 


For that reason we will investigate the post-Newtonian expansion in if. Interestingly 
the convergence for the first nine coefficients at 7 = 1 is faster in if then in e. The 
extreme Kerr-Newman black hole obeys the parameter relations 


MO (2 - if 2 ) 

\A 


1, 


(48) 
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J 2 


and 


(1 - if 2 ) M 4 
a\U) = 


= 1 , 




(49) 


(50) 


2 - V’ 2 ' 

If there is a transition to the extreme Kerr-Newman black hole, then certain 
functions R( if) must go to zero for 7 —>■ 1 and if E [0,1]. In that case with 

(51) 


lim 11 ( 7 ) = 0 , 
7— >1 


the inequality 


lim S'( 7 , if) 

7—>1 





7—U(7) 


or 


lim S 1 

7—>1 


(7,^) 



lim y (7 \ 


701 R( 7 , if) 


<00, ^ e [0,1] 


> 0 , if e [ 0 , 1 ] 


(52) 


(53) 


hold for a particular i'(y). It remains to show that S{^ : if) or S 77 ,if) converge at 


7 = 1. For that we can use the convergence criterion (41) for the majority of cases. 


In the majority of cases, we will present the results in a split figure. In part (a) we 
present the value of a quantity for 7 = 1 using increasing expansion orders N to show 
the convergence behaviour and sometimes using a Pade approximation in y assuming 
a better convergence then the power series expansion. I 11 part (b) we present the values 


of the convergence criterion (41) for increasing orders n. The critical value 1 is marked 
by a dotted line (red). 


Functions on the disc: First of all we want to show that the condition (43) is fulfilled. 
This follows if f'(g* < 1, (* = 0) = 0. For that purpose we investigate 

Sn = [fie* < 1, C* = 0)] 1/2 (1 - 7) _1/4 = 1 - (j - ) 7 + 0 (7 2 ) (54) 

In Figure [I] (a) we see that S'd stays finite on the disc for 7 —* 1. Figure [I] (b) shows 


that the convergence criterion (41) is fulfilled. For the sake of clarity, we plotted only 
Sb( 7 =l, JV=10) and K 9 {S D }. 

The minimum of K 9 {»Sd} is at if = 1, g* ~ 0.8. The convergence criterion €3 is 
also fulfilled by KyjjSn} and K 8 {Sb}. The shapes of the plots of KsIS'd} up to K 8 {Sb} 
are similar to 1K 9 {S'd} with lower values. So IK^IS'd} is increasing with higher values 
of n. We plotted the ECD case (if = 1) in Figure |2j were the convergence is at worst. 
The values for Sj)(if = 1,7 = 1,1V) are decreasing with increasing orders N. The Pade 
approximation gives a much better result. The convergence criterion is fulfilled. Note 
that in the ECD case it is sufficient for a black hole limit that go > 0 for 7 —y 1. 

We have the analytic value Sd(p*=0) = Sn(if= 0 ) = (1 — 7 ) 3//4 , and this coincides 
with Figure [I] (a) and the values Sb(p*=0,7=1, N— 10 ) = Sn(if= 0 , 7=1,1V=10) = 
0.048568.... 


Altogether, we find f(g* < 1,C* = 0) = 0 and with (34) that (43) is fulfilled. 
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Figure 1: The quantity [f(g* < 1,£* = O )] 1 ^ 2 (1 — 7 ) 1/4 in the limit 7 —» 1 for the 
expansion orders N = 10 (a) and the convergence criterion for n = 9 (b). 



(a) (b) 

Figure 2: The quantity < 1, £* = O )] 1 ^ 2 (1 — 7) _1/4 for the ECD case if) — 1 in the 
limit 7 —>■ 1 for increasing expansion orders N and for a Pade approximation P[10,10] 
in (dashed dotted line) (a) and the convergence criterion for increasing n (b). 
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Disc radius: In order to show, that the disc radius go vanishes for 7 —> 1, we will 
investigate the ratio 


M 


h'M = — 1 (!- 7 ) = tti 7 
Qo on 

In Figure |3] (a) we can see that the inequality 


1 + ( 15 “ tA* ) 7 + ° O) 


(55) 


iim nj > 0, 

7 —>1 


e [ 0 , 1 ] 


(56) 


of type (53) is fulhlled. Figure [3] (b) shows, that the K n from (41) exceed one for large 
values of n. 

The analytic value in the uncharged case is 


^(7 = 1,^ = 0) = 0.32863..., 

while the Pade approximation (blue dashed dotted line in figure [3] (a)) gives 
P[6,12 ]{/^( 7 = l,^ = 0,iV=9)} = 0.32885.... 


We also expanded hm in 7 and plotted the inverse Pade approximation 
(P[6,12]{/i M (7 = l)}) -1 as a green dashed dotted line in figure [ 3 ] (a). To compare 
the values, we have: 

(P[ 6,12]{/i m ( 7 = 1,$ = 0,iV=9)}) 1 = 0.32872..., 

(P[6,12]{/i M ( 7 = l,if = l,i\T=9)})- 1 = 0.19675..., 

P[6,12]W(7 = 1,-0 = 1, N= 9)} = 0.20117.... 


All in all we find that in the limit 7 —y 1 the disc radius po tends to zero. 
Parameter relations from the extreme Kerr-Newman black hole: Now we check the 


parameter relations (48) and (49). Note that Q and J have the global pre-factor 


y/l — fj 2 . At first we will investigate the quantity 

Mfi(2-?/; 2 )l. i 1 / 8 4 (2 \ 2 , 

(!- 7 ) =l+(l-^ + ^0T +0(7”). 




(57) 


77(1 — V 2 ) 

In Figure |4] (a) we can see that Sq stays finite at 7 = 1. The singular behaviour of the 
around if = 0.5 (b) can be explained by the zeros of the coefficient functions c n (c). 


Although (41) is not valid in a small interval in fj, the convergence of Sq is extremely 


good there, which can be seen in (a). 
It turns out that the quantity 


S.j = 


1 

2 + 


1 - 


J° 


t/tIT 



(1-7) 


-1 


(1-7) 


-1 


37T 

1 + 

f 39 5 

1 

107 

^20 7T 

14 


r 7 + 0(7 2 ) 


(58) 

(59) 
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(a) (b) 


Figure 3 : The ratio M*( 1 — 7) in the limit 7 —> 1 for increasing expansion orders N and 
for two Pade approximations P[ 6, 12 ] in 77 (dashed dotted lines, blue/green) (a) and 
the convergence criterion for increasing n (b). 


converges and stays finite in the limit 7 —> 1 (see Figure [ 5 ]). This is a check for ( 49 ) and 
moreover it means that J° can be written as 


J° = \A-V> 2 [1 + (C 2 - 3/8) (1 - 7 ) 2 + O (1 - 7) 3 ] , 


with lim Sj = —C2. We will use this result in section 

7 —>1 


0 


If ( 48 ) and rt 49 ft hold for the disc at 7 = 1 , it follows with ( 35 ) that 


lim of = 

7—>•! 




2 


In accordance with ( 50 ) we identify the constant in ( 43 ) to be 


lim a(g* < 1 , C* = 0) = if (2 - ^ 2 ) 


2 \ —1 


( 60 ) 


( 61 ) 


( 62 ) 


Conclusion: Within the high accuracy of these calculations we can state that the 
solution describing rigidly rotating discs of charged dust has a transition to an extreme 
Kerr-Newman black hole. The multipole moments also fit in this result, as implied in 
the next section. 


5. Leading order behaviour close to the black hole limit 

Interestingly, it turns out that the power series expansion at 7 = 0 could be used to 
make a statement about the leading order behaviour for the power series expansion 
at 7 = 1 . The uncharged case is discussed analytically in [ 20 ]. For this purpose we 
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(a) (b) 



(c) 


Figure 4: The quantity Sq in the limit 7 —>■ 1 for increasing expansion orders N and for 
a Pade approximation P[ 4,16] in y'y (dashed dotted line) (a), the convergence criterion 
for increasing n (b) and the coefficient functions |c n | (c). 


investigated the coefficients from (24) in two normalized and dimensionless forms: 

m r 


= 


m 


n+1 


and e° = 


m. 


n+l 


with mo = M, mi = iJ, 


(63) 


* m n / 1 , 2\ n / 2 


and 


e* = 




2\ n / 2 


with k n = 


(64) 


k n f' k r 

Herein the k n correspond exactly to the coefficients m n and accordingly e n /"0 of a 
Kerr-Newman solution with mass M, angular momentum J and charge Q. Thus 
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Figure 5: The quantity 1/2 + 1 — J°/ -y /7 (1 — ip 2 ) (1 — 7 ) 1 (1 — 7 ) 1 in the limit 
7 —>• 1 for increasing expansion orders N and for a Pade approximation P[ 3,15] in y /7 
(dashed dotted line) (a) and the convergence criterion for increasing n (b). 




(a) (b) 


Figure 6 : The quantity yjl — fj 2 — e^ (1 — 7 ) 2 in the limit 7 —>■ 1 for increasing 
expansion orders N and for a Pade approximation P[13, 5] in ^ (dashed dotted line) 
(a) and the convergence criterion for increasing n (b). 


m n = e* = (1 — for black holes, and they are equal to m° n and accordingly 

e 0 n /f> in the extreme case. 
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In the black hole limit we checked 


dm” 


d7 


7=1 


de° 

d7 


= 0 , 


dm* 


7=1 


d7 


7=1 


de* 

d7 


= 0 


(65) 


7=1 


up to n — 11 , while the second derivatives of the m* at 7 = 1 are only zero in the 
uncharged case. We already showed that (65) holds for m°. It follows with (60) that 


dm” 


dy 


= 0 yy 


dm* 


7=1 


dy 


de° 

= 0, —^ 
7=1 d 7 


= 0 yy 


7=1 


de* 

dy 


= 0 


( 66 ) 


7=1 


for any fixed value of n. This is important, because some quantities have a better 
(faster) convergence then others. 

We present the following two important quantities as examples: The magnetic 
dipole moment e\ and the gravitational quadrupole moment rn ? 2 , namely we investigate 


S ei = 


and 


\/l -V ~ e\ (l-y) 2 = 


1 + 


32 1 

- ib' 

35 5 r 


S m2 = [i 2 (1 - ip 2 ) - ml] (l-y) 2 = 


97T 2 

8 O 7 2 


1 + 


23 18 

- ip z 

35 35 ^ 


1 + 0 ( 7 2 ) 


7 + 0 ( 7 2 ) 


(67) 


( 68 ) 


ei and S m2 stay finite in the limit 7 —» 1 . The 


In Figure [ 6 ] (a) and [T] (a) we see that S, 
convergence criterion for S ei is fulhlled in the regular parameter space and for S m2 up 
to \ip\ < 0.99 at n — 8 (see Figure [7] (b)). We hnd a similar convergence behaviour for 
the electric quadrupole moment e 2 , while for higher multipole moments, up to n = 11 , 
the convergence criterion is fulhlled in the entire regular parameter space. 

For the following discussion we assume that (65) holds for all n > 0. If this is the 
case the Ernst potentials £ and $ can be expanded at 7 = 1 in the following way: The 
series starts with the Ernst potentials £kn an d ‘Fkn of the Kerr-Newman solution with 
mass M , charge Q and angular momentum J and the first order term in (1 — 7 ) is zero: 


S(M, J, Q- 0 = T kn (M, J; C) + O [(1 - y) 2 ] , 

(69) 

, (v ,. n M(C-M)-i./ 

(70) 

$(M, J, Q; 0 = $kn(M, J, Q- () + O [(1 - y) 2 ] , 

(71) 

#kn(M,J,Q,C)- mk + m) _. / . 

(72) 


Comparing (60) with the power series expansion from J° at 7 = 0, we can state that 

T e [0,1], e G [-1,1]. (73) 


It is equal to one ify = lor^ = ±l. This means that £kn an d ^kn ar e given by 
the axis potential of a “hyperextreme” Kerr-Newman solution for 7 < 1 and |^| < 1. 
The Kerr-Newman solution already contains the first multipole moments M, J and Q. 
Therefore, in the far held the second order term in (1 — 7 ) for £ starts with order 
whereas the term for starts with order <C~ 2 . 
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(a) (b) 

Figure 7: The quantity [i 2 (1 — if 2 ) — rn°] (1 — y )~ 2 in the limit 7 —>■ 1 for increasing 
expansion orders N and for a Pade approximation P[15, 3] in ^7 (dashed dotted line) 
(a) and the convergence criterion for increasing n (b). 


All together we found strong evidence that the disc solution near the ultra- 
relativistic limit at 7 = 1 , and not to close to the disc itself, could be approximated very 
well by a “hyperextreme” Kerr-Newman solution with the same gravitational mass, 
angular momentum and charge. 
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Appendix: Coefficient functions 


General structure of an approximated function (see (40)): 


F('y,if,N) = F N l(7,V0 


N 


1 + ^ C nidf)l r 


n =1 


(•i) 


• Charge parameter: if — Q/M — e y 1 — E^j . For the relative binding energy 

see Appendix C. in [5]. 

• Funtion on the disc f'(g* < 1, £* = 0): 

/nl = 1 ) 
ci = —2 + if 2 g 2 , 


(• 2 ) 

(.3) 
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c 2 = 1 - (^ 2 - -/) ^ 

V2711 16 \ 2 /1027 16 

Cs “ V 3200 ~ 9tO ) g + V 2560 OO 


3072 9tt 2 


* 


f 29713 8 \ 2 

\ 22400 ~ Tk 2 ) 6 + 

- 4) W 

33600 9n 2 n ^ 


1633 16 

3840 ~~ OO 


631 16 

4608 ~~ 9tO 


-Ah> 6 ^ 


Gravitational mass M *: 

M* = — 

NL 3vr’ 

9 1 /2 

10 15 y 


_ 11 64 

C 2 ~Y~9tO 


58333 8 \ , 2 

67200 “ + 


1429 8 W 

33600 9tt 2 J ’ ’ 


2521 1568 

1260 ~~ 135tt 2 


393917683 166013 

344064000 ~~ 15120tt 2 


301426021 3175 

1548288000 ~~ 3024tt 2 


22109 56 

756000 ~~ 135tt 2 


6242 23792 | 4096 

3465 ~~ 945tt 2 + m 4 


1725516795817609 1658106057221 1024 

2790386565120000 ~~ 40874803200tt 2 + ~3tO 


4757113193121521 | 2400073720591 128 

9766352977920000 + 122624409600vr 2 ~ TOO 
91538340739 2245701613 128 \ /6 

1300561920000 ~~ 479001600vr 2 + 3 ^) ^ 
11617705117 1907 64 \ 

447068160000 3780vr 2 + TOO ) ^ ' 


Angular velocity 0° = QM: 

r'vo 4 7 3 / 2 r -— 


^nl — v4 l/’ 2 ; 

3 3 /2 

Ci =-- xb , 

20 10 V ’ 


1249 32 

1120 ~ 3tO ' 
1219 280 


9433 12 \ , 2 

8960 ~~ tt 2 j ^ + 


14447 4 

67200 ~ 3t^ 




1219 280 / 18872669 5681 \ /2 / 27689221 

1152 ~~ m 2 ~~ \ 294912000 ~ 4320vr 2 ) >P ~ V 30965761 
/ 988499 22 \ /6 

+ V 6048000 ’ 

4861867 29993 6656 

7096320 ~~ 945tt 2 + 200 

( 393976735547767 296409922049 1664 

+ V 1860257710080000 + 5449973760vr 2 ~~ ~3tO 


2768922671 8323 

3096576000 ~ 864vr 2 
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117736707196657 275650267021 10088\ /4 

6510901985280000 + 7431782400vr 2 ~~ 27 tt 4 ) U> 
7449051923917 4791322541 208 \ /6 

9537454080000 ~~ 319334400vr 2 + 3 ^) ^ 
27346973 53287 104 \ /8 

197120000 ~~ 30240vr 2 + TOO / ^ ' 


Angular momentum J° 
J° = 

NL 10^7 v ’ 

1 1 / 2 
Ci =- w , 

20 14 7 

77363 416 

C2 “ ~50400 + TOO ^ 


77363 416 7162173 308 \ , 2 7 3211 

_ 50 4 00 + TOO + V134400 ~~ TOO ) ^ V 86400 

361961 488 7 4379136107 429269 \ ; 

~336000 + 400 + V10597171200 ~~ 110880vr 2 ) ^ 
7145460164961 683051 \ ;4 

+ V 238436352000 ~~ 110880vr 2 ) ^ 

7 11652667 94 \ /6 

~ V 465696000 ~~ 315tt 2 ) ^ ’ 

5382534221 | 10753 14848 

3228825600 + 4725tt 2 ~~ 81vr 4 

7619674800505859183 57472167761 


86400 9tt 2 / V ’ 


28288 


169283451617280000 70849658880vr 2 81 tt 4 7 v 

135571135157430463 3757347355969 15464 \ 4 

84641725808640000 + 1062744883200tt 2 ~~ 81vr 4 ) ^ 
2236889283151 8830544257 2096\ /6 

5904138240000 ~ 1383782400vr 2 + 800) ^ 

206530297 135943 8 \ /8 

10567065600 453600vr 2 + 9tO ) ^ ' 


Magnetic dipole moment D°: 
37T , r -— 


29 9, 2 

C 1 - ™ ™ ^ > 


28 70 T 7 

13099 656 

Co —-“f" -7 


c 3 = - 


13099 656 71637029 428 \ , 2 7 9 4 \ , 4 

~ 5600 + TOO + V1209600 ~~ m 2 ) ^ + \ 86400 + ^0) ^ ’ 
14126969 56044 779899457229 5901503 \ /2 

~ 5174400 + 2079 ^0 + V 43352064000 ~~ 332640vr 2 ) ^ 

7143489392481 2100353 \ /A 7 326299 22 

+ V 238436352000 ~~ 332640vr 2 / ^ \ 42336000 ~~ 105 tt 2 

104788527749 | 20606291 22528 

80720640000 + 1351350vr 2 ~~ 81 tt 4 


W-24) 
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/2217555015040271693 2568918307691 


40768 


V 507850354851840000 ' 354248294400vr 2 81 tt 4 

/ 45857036770964723 1399751492287 20744\ 

v 16928345161728000 ~~ 1062744883200vr 2 ~~ 81 tt 4 J 
f 16546156571327 25927117379 2576\ /6 

V 53137244160000 ~~ 4151347200vr 2 + 817T 4 J ^ 

f 7617384713 29677 8 \ /8 

151200vr 2 + 9^ ) ' lp ' 


v968647680000 
Gravitational quadrupole moment mf. 

O _ 9 

m 2NL ~ 


807r 2 7 2 ’ 


C 2 = 


C 3 = 


47 

_ 35 ~ 
4267 

±v, 

35 ^ ’ 

464 / 

1 1 1 

4 3667 

512 N 

6300 

+ 27tt 2 + 1 

v 2240 

27vr 2 , 

12371 

256 

/ 9772341947 

40425 

105tt 2 

V39739392000 


^ 2 - 


^ 22237 


16 


—- 1 

151200 9tt 2 1 ’ 


50123 \ 


( 81384847 . 113 \ , 4 


r- 


c 4 = 


+ V 3973939200 + 5544tt 2 J 
11295719 128 4096 

6006000 ~~ 457T 2 TOO 

(5428030087901639 7641678347 


/ 3465793 
\58212000 


1024\ 


V1692834516172800 2683699200vr 2 3 tt 4 

62°8\ 
27vr 4 ) V 




I' 69407086543177 48426129893 

+ 5904138240^2 


v46506442752000 
( 3595050520793 1093286657 128 \ 

V 30996725760000 + 345945600tt 2 + 3 ) 
( 13554494567 


V 322882560000 18900vr 2 

Electric quadrupole moment e\\ 

9 


12287 64 \ f 

+ 27^4 ^ 


e 2NL — 


807r 2 7 : 




59 6 o 

c i =- ib , 

35 35 V ’ 

673 464 

c 2 = 


2100 27vr 2 

9575749 10112 


/78401 
{60480 


512 ) 


c 3 = 


9702000 


27tt 2 
/23595907079 
l 17031168000 


f 2741 
~ V16800 
17119 \ 


16 \ /4 


r 


c 4 = 


r 


64 \ 
105t r 2 J 




11557T 2 V 17031168000 1320vr 2 / 

/1433090333 3159 \ ; 4 / 10037 32 \ /6 

+ V 2838528000 ~~ 616tt 2 / ^ \ 117600 ~~ 35tO ) ^ ’ 

16364265889 1388 4096 

10090080000 ~~ 5005tt 2 ~ TOO 


(.25) 

(.26) 

(.27) 

(.28) 

(.29) 


(.30) 

(.31) 

(.32) 

(.33) 

(.34) 
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265686220800?r 2 
436389127799 


3tt 4 J ^ 
6208 ^ r 


(32617425522286343 1886067001921 

v11542053519360000 + 

/25607911281592511 

\21160431452160000 + 37955174400?r 2 ~~ 27 ^) 
f 208667707711 467429993 128 \ /6 

V 3444080640000 ~~ 94348800vr 2 + 3 ^) ^ 
f 1019676659 2249 64 \ 


V 17611776000 2700vr 2 


+ 


27tt 4 J 


r- 


(.35) 
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